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We present an exact field theoretical representation of an ionic solution made of charged hard 
spheres. The action of the field theory is obtained by performing a Hubbard-Stratonovich trans- 
form of the configurational Boltzmann factor. It is shown that the Stillinger-Lovett sum rules are 
satisfied if and only if all the field correlation functions are short range functions. The mean field, 
Gaussian and two-loops approximations of the theory are derived and discussed. The mean field 
approximation for the free energy constitutes an exact lower bound for the exact free energy, while 
the mean field pressure is an exact upper bound. The one-loop order approximation is shown to be 
identical with the random phase approximation of the theory of liquids. Finally, at the two-loop 
order and in the pecular case of the restricted primitive model, one recovers results obtained in the 
framework of the mode expansion theory. 
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I. INTRODUCTION 



Various ionic systems including electrolyte solutions, molten salts, and colloids can be studied with a good approx- 
imation in the framework of the so-called primitive model (PM) which consists in a mixture of M species of charged 
hard spheres (HS) which differ by their respective charges and (or) diameters. 1 Of special interest is the restricted 
primitive model (RPM) where M = 2, the hard spheres have all the same diameter, and the cations and anions bear 
opposite charges ±q. In many instances, we shall also consider the special primitive model (SPM), where the number 
M of species as well as the charges are arbitrary but all the ions have the same diameter a. 

In the two first parts of this work, published some years ago and hereafter referred to as I and II, we have established 
an exact field theoretical representation of the RPM. 2 ' 3 The action of this field theory, which is obtained by applying 
the Kac-Siegert-Stratonovich-Hubbard-Edwards (KSSHE) 4-9 transform to the Coulomb potential, looks like the sine- 
Gordon action to which it reduces in the limit of point-like ions, 10-12 hence the slightly abusive title of this series of 
papers. Nowadays we prefer the acronym KSSHE to christen the action. The extended sine-Gordon action derived in 
paper I for the RPM is obtained here for a general PM. The regularization of the Coulomb potential which is required 
to define properly the KSSHE transform is obtained by a smearing of the charge over the HS volume. A more general 
treatment where a part of the Coulomb interaction is incorporated in the reference system is discussed in the review 
of Brydges and Martin. 11 

The developments of refs. I and II are based on a cumulant expansion of the grand partition function reorganized in 
ascending powers of either the fugacity or the inverse temperature. In this way one can obtain the exact low fugacity 
and high temperature expansions of the pressure and the free energy of the RPM. Of course the expressions obtained 
in that manner are already known from the theory of liquids and were derived years ago in the framework of Mayer 
graph expansions. 13,14 In the present paper we proceed differently. After having obtained the KSSHE action for the 
general PM (see sec. II), we reorganize the cumulant expansion by grouping some classes of Feynman diagrams. The 
resulting loopwise expansion is explicitly computed up to order two in the number of loops. In the case of the RPM, 
the two-loop order free energy turns out to coincide with an expression derived more than thirty years ago by Chandler 
and Andersen 15 in the framework of the so-called mode expansion theory. 16 ' 17 Reorganizing the loop-expansion in 
ascending powers of the inverse temperature gives back the high-temperature expansions of paper II and ref. 
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Our paper is organized in the following way. In next section II we show how to construct a well-defined KSSHE 
transform by regularizing the Coulomb potential at short distances by means of a smearing of the charges inside the 
volumes of the HS. In section III we establish the general relations between the charge correlation functions and the 
correlations of the KSSHE field. From the known asymptotic behavior of the former one can deduce that of the 
latter. The conclusion, which is detailed at length in section IV, is that the n-body correlations of the KSSHE field 
are short-ranged functions; stated otherwise, the KSSHE field is a non-critical field. The so-called Stillinger-Lovett 
sum rules, both for the homogeneous 18 and the inhomogeneous fluid 19 emerge as a consequence of this behavior. In 
section V the mean field (MF) level of the theory is studied in detail. The MF free energy (3Amf is shown to be a 
strictly convex functional of the M partial densities and to constitute a rigorous lower bound of the exact free energy. 
The former property excludes a fluid-fluid transition at the MF level while the latter serves to define an optimized MF 
free energy by maximizing (5Amf with respect to the variations of the smearing functions. An explicit expression of 
the optimized (3Amf is obtained in the case of an homogeneous fluid. From the MF solution for the inhomogeneous 
system we also deduce the expressions of the n-body correlation and vertex functions of the homogeneous system in 
the Gaussian approximation. 20 This Gaussian approximation is discussed in section VI and shown to be equivalent to 
the random phase approximation (RPA) of the theory of liquids. 1,21 Finally, a two-loop order calculation is performed 
in section VII. The resulting expression for the the free energy of the RPM is shown to be identical with that obtained 
by Chandler and Andersen 15 in the framework of the first version of the mode expansion theory. Conclusions are 
drawn in section VIII. 



II. THE KSSHE TRANSFORM 



A. The model 



We shall consider only the three dimensional (3D) version of the (PM), i.e. a mixture of M species of charged hard 
spheres. 1 The ions of the species a (a = 1, ... , M) are characterized by their diameter a a and their electric charge 
q a . The molecular structure of the solvent is ignored and it is treated as a continuum, the dielectric constant of which 
has been absorbed in the definition of the charges q a . The solution is made of both positive and negative ions so 
that the electroneutrality in the bulk can be satisfied without adding any unphysical neutralizing background to the 
system. The particles occupy a domain f2 C M 3 of volume fl of the ordinary space with free boundary conditions. 
Only configurations w = (N\; f\ , . . . , f 1 Ni \. . . \N a ; rf , . . . , fjy \Nm; rf 1 , . . . , rj\f M ) (rf € ft ) without overlaps of the 

spheres - i.e. such that \\rf a — rf | > (o~ a + (7p)/2 - do contribute to the partition or grand partition functions. In such 
a configuration, the charge q a of each ion can be smeared out inside its volume according to a spherically symmetric 
distribution q a T a (r) without altering the configurational energy as a consequence of Gauss theorem. The distribution 
r a (r) is a priori arbitrary, provided it satisfies the following properties : 



/ 



T a {r) = if r > a a = a a /2 , (2.1a) 
d 3 rr a (r) = 1 . (2.1b) 



The electrostatic interaction energy of two charge distributions r Q and the centers of which are located at the 
points r\ and r 2 of respectively will be noted w a ^(l, 2). It reads as 

w a , p (i,2) = J d\ v J d\ 2 , r^n - M)v c (\\n> ~ M)T0(\\f 2 > - f 2 \\) , 

eet q (1,1')« c (1',2')t /3 (2',2) , (2.2) 

where v c (r) = l/r is the Coulomb potential. Note that in this paper, summation over repeated, either discrete or 
continuous indices will always be meant (except if explicitly stated otherwise). As a consequence of eqs. (2.1) and of 
Gauss theorem w a .p(l, 2) = l/ri2 for r\ 2 > a a + cfp. Note that the Fourier transform w a ,(3(k) of the interaction takes 
the simple form 

which diverges for k — > as Air/k 2 since r Q (0) = 1, as follows from eq. (2.1b). Finally we shall denote by 

v a ,p(l, 2) = q a qp w a ,p{l, 2) (2.4) 
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the pair interaction of two ions. 

The electrostatic potential energy of the configuration ui times the inverse temperature 13=1/ kT can be written as 

fiUaiu) = | p c (l) « c (l, 2) pc(2) - N a u s a , (2.5) 

where /3c (1) is the microscopic charge density in the configuration uj at the point f\ and is the self-energy of the 
charge distribution q a r a (r). In general, for a sufficiently regular distribution r a (r), the self-energy 

= ^ »a,«(0) = ^ jd-k p a (kf , (2.6) 

where d v /c = ci 3 fc/(27r) 3 , is a well-defined positive and finite quantity. Of course diverges for point-like charges which 
makes the KSSHE transform, to be introduced in next section, an ill-defined object in that case. The microscopic 
smeared charge density pc(r) which enters eq. (2.5) reads 

Pc(l) = «aT Q (l,l , )?a(l'). (2-7) 

where 

N a 

Pa(l)= E * 3 (*-*?J ( 2 - 8 ) 

*a = l 

is the microscopic number density of the species a at the point f\. 

It will prove convenient to make use of Dirac's notations for matrix elements and scalar products and to rewrite 
the energy (2.5) as 

ipc(lK(l,2)j3c(2) = i (pc\v c \pc) , 

= (pc\v) , (2.9) 

where V(\) = pc(l') w c(l', 1) denotes the microscopic electric potential at the point f*i in the configuration u. Of 
course V is solution of the 3D Poisson equation, i.e. 

Ai V(l) = -4wp c (l) . (2.10) 



B. The KSSHE transform of the Boltzmann factor 

The Boltzmann factor in the configuration u> is equal to 

exp {-0U H) = exp {-f3U HS (w)) x cxp (-(3U el (w)) , (2.11) 

where Uhs (w) denotes the contribution of the hard cores to the configurational energy. We perform now a KSSHE 
transform in order to rewrite eq. (2.11) as 2-12 

exp (-0U M) = cxp (-pU HS M) cxp (N a v s a ) (cxp (i/3 1 / 2 (p c \<fi)) ^ , (2.12) 

where the brackets (. . .) denote Gaussian averages over the real scalar field <p(r), i.e. 

<■■■>„„ =N- 1 jv V ...e^^- l -( V \v- l \^ , 

K e = J ^cxp (-± (vlv- 1 ^ , (2.13) 

where 
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« c - 1 (l,2) = -i-AiJ(l,2) (2.14) 
is the inverse of the positive operator v c (l, 2). Therefore one has, after an integration by parts 

M Vc = fvip cxp (-1- d 3 r (V^) 2 ) (2.15) 

The functional integrals which enter eqs (2.13) and (2.15) can be given a precise meaning when grounded perfect 
conductor boundary conditions (BC) are adopted; periodic BC's also work if only neutral configurations are considered, 
we refer the reader to the literature for more details. 2,11,20,22 ' 23 It will be convenient to write 

(fi C \<p) = {pM , (2.16) 

where the smeared field 4> a is defined as 

Ml)^?/ 2 q a T a (l,l')tp(l') . (2.17) 

The field i<j> a (l) may thus be seen as an external one-body potential acting on the particles of the species a; indeed, 
one can rewrite the Boltmann factor (2.12) under the form 



exp (-(3U (u)) = cxp {-/3U H s M) cxp (N a v% ) (cxp {i{p a \^ a )))v c , 

/ M N a \ 

= exp {-/3Uhs M) cxp (N a i/f ) (exp ]T ]T i<t> a (r?J ) Ve . (2.18) 



\a— 1 i a — 1 



C. The Physical meaning of the auxiliary field 

In a given configurational ui let us define an action 

h[<p] = \{<p\Vc 1 \<P)-i (Mpa) , (2.19) 

and a partition function 

z{u>) = (exp{-h[<p])) Vc . (2.20) 

(Henceforth we shall specify the arguments of functionals by means of brackets, the variables of ordinary functions 
being enclosed as usual by parenthesis.) The saddle point of the functional h[ip] is obtained by solving the equation 



Sh 



_ = 0, (2.21) 



5ip(r) 

which can be recast under the form of the Poisson equation 

Ai5>(l) - -Am/S^pcil) , (2.22) 

the solution of which is of course 

^(1) = z/3 1/2 t>(l) . (2.23) 

Therefore, at the saddle point, the field Tp can be identified with the microscopic electric potential in the configuration 
u, up to an imaginary multiplicative constant. Moreover, it is easy to show that the value of h[ip] is nothing but the 
energy of the configuration, i.e. 

%] = f (PcK\pc) ■ (2.24) 

Let us make now the change of variables tp = Tp + Sip where Sip is a real scalar field. It follows from the stationarity 
condition (2.21) that 



4 



h[<p] = h[ip] + ±(5<p\v- 1 \6<p) , 



(2.25) 



which confirms that h\0\ is indeed a minimum of the functional h[ip] since v~ is a positive operator. For a more 
complicated Hamiltonian than h[<p], the approximation consisting in truncating its functional Taylor expansion about 
the saddle point at the second order level is called the Gaussian approximation. 20 This approximation is obviously 
exact for z(ui) because h[ip] is a quadratic form. A direct calculation indeed confirms that 

/ \ ( ,r- n J^exp^i^l^- 1 !^)) 
z(uj) = exp(— h[ip\) 
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= exp(-|<?cH?c>) • (2-26) 



D. The KSSHE transform of the grand partition function 

Henceforward we shall work in the grand canonical (GC) ensemble. We denote by fj, a the chemical potential of the 
species a and by ip a (f) the external potential with which the particles of the species a interact eventually. According 
to a terminology due to J. Percus, 24 we shall define the local chemical potential v a (f) as (3(fi a — ipa{r))- With these 
notations, the GC partition function of the system takes the form 

= E am • • • E ATI / d ^ • • ■ d ^N M exp(-/3C/(o;)) jj <=x P {v a (rfj) . (2.27) 

Grand canonical averages of dynamic variables A(u)) will be noted (A(u))) GC . Inserting the expression (2.18) of the 
Boltmann factor in eq. (2.27) one obtains the KSSHE representation of S 

- (E HS [{V a + i<j> a }]) Va , (2.28) 

where V a = v a + and Ehs [{v a + i4>a}] is the GC partition function of a mixture of bare hard spheres in the 
presence of the local chemical potentials V a + i<fi a . The above result generalizes to the case of the PM the result 
obtained in paper I for the restricted primitive model. It is also possible to incorporate a part of the Coulomb 
interaction in the reference potential which yields a more general expression than cq. (2.28) as detailed in the review 
of Brydges and Martin (cf eq. (2.29) of ref 11 ). However, in the liquid domain, the thermodynamics and correlations 
of this reference system are, by contrast with those of the HS fluid, little known in general. Relations similar to eq. 
(2.28) have also been obtained and discussed for neutral fluids. 5,8 ' 25 ' 26 

To make some contact with statistical field theory we introduce the effective Hamiltonian (or action) 

W<P\ = ^K 1 !^) - logE HS [{V a +icj> a }} , (2.29) 

which allows us to recast 5 under the form 

E[{u a }] = AT' 1 J exp(-WM) . (2.30) 

It will be important in the sequel to distinguish carefully, besides the GC averages < . . . >gc, between two types 
of statistical field averages : the already defined < . . . > Vc and the < . . . >u that we define as 

< m >ns fJgjM -"MMM . (2 . 31) 

J Vtp exp{-H[<p\) 

With these definitions in mind one notes that for an arbitrary functional A[f] one has the relation 

.4 r -I i A [<P] E HS [{"a + i^a}]) Vc 

< A[ip] > n = — p= ,j- - 7-7 m ■ (2.32) 



5 



III. CORRELATION FUNCTIONS 



A. Zoology 

The ordinary and truncated (or connected) density correlation functions will be denned in this paper as 27 ' 28 

G^.a B [K}](l,...,n) = (np ai (i)) , 

\l=l / GC 

= *\<u 6nE[{Ua}] 



Our notation emphasizes the fact that the G$... an (truncated or not) are functionals of the local chemical potentials 
v a (r) and functions of the coordinates (1, . . . , n) = (f\, . . . , r n ). Note however that, in the remainder of the paper, 

(n) 

we shall frequently omit to quote the functional dependence of G ai ... a „ upon the v a when no ambiguity is possible. 
In standard textbooks of liquid theory 1 the n-body correlation functions are more frequently defined as functional 
derivatives of S or logS with respect to the activities z a — exp(j/ Q ) rather than with respect to the local chemical 
potentials. This yields differences involving delta functions. For instance for n — 2 and for a homogeneous system 
one has 

^a/H 1 ' 2 ) = PaPp9a(i(ri2) + PaS a ,/3 5(1, 2) , 

Gap T (!> 2 ) = P a Pph a p(r 12 ) + p a S a . p 5(1, 2) , (3.2) 

where p a is the equilibrium number density of the species a and g a p(r) the usual pair distribution function; finally 
h a p = g a p - 1- 

The charge correlations will play an important role in subsequent developments. They are defined as 

GP(l,...,n) = /f[pc(i)) ■ (3-3) 

\l=l / GC 

It follows from the definition (2.7) of the smeared density of charge pc that eq. (3.3) can be rewritten alternatively 
G%\l, ...,n) = q ai ...q an r ai (1, 1') . . . r Q „ (n, n') (l', . . . , n') . (3.4) 

Clearly the operator 

e(l) = z/3 1 / 2 g ct r Q (l,l') -4-77 (3.5) 
is the generator of the charge correlations for we have clearly 

%n/3 n/2 G (n) (1; _ . _ ; n ) = S -l (l) . . . 0(„)S . (3.6) 

The truncated charge correlations can thus be defined according to 

inpn/2 Q{n) T (1) ^ ^ = _ _ _ Q(n) Jog g _ (g ?) 

On the one hand 

GP T (l, . . . , n) = q ai . . . q an r ai (1,1')... r a „ (n, n') G^ T an (l', ...,ri), (3.8) 
and, in the other hand 22 ' 27 ' 28 

G%> T (l, . . . , n) = GP T (l,...,n)-Y J Y[G { ™ )T ^,...^ m ), (3.9) 



G 



where the sum of products is carried out over all possible partitions of the set (1, . . . ,n) into subsets of cardinality 
m < n. The functions G^ ' (resp. G^ ' T ) for different values of n are not independent; they are related by a hierarchy 
of equations most conveniently written with the help of the operator defined at eq. (3.5). The hierarchies for the 

G { c ] and the G^ 7 are derived in appendix A. 

In the field theoretical representation of the PM the field correlation functions play a key role. They are defined as 

GW(l,...,n) = <¥>(!) •■■¥>(«)>« . (3-lOa) 
GW T (l,..,n) = GW T (l,..,n)-^ J] r (ii, . . . ,i m ) • (3.10b) 

Of course the G^\ as the charge correlation functions, are functionals of the local chemical potentials. The hierarchies 
for the G^ and the G^ T are derived in appendix A. 

B. Relations between the charge and field correlation functions 

1. The density and charge correlation functions as statistical field averages 

It follows from the definition (3.1) of G^... a „ and from the KSSHE representation (2.28) of the grand partition 
function that we have 

os..„. Mia «) » P (->-„>) i^'^i' . 

= H-'(H„ S C« „,...„„ [{F„+i*,}](l,...,»))^ , (3.11) 

where G^g a a [{V a + i<f> a }] (1, . . . , n) denotes the density correlation function of the reference HS fluid in the 
presence of the local chemical potentials {V a + i(f> a }. Thence, making use of eq. (2.32) 

Git. an [K}](1, ■■■,«) = (g ( hI „,..«, [{^ + i0a}] (1, • • • , • (3-12) 

Eq. (3.12), which extends to ionic mixtures a relation that we derived elsewhere for simple non-charged fluids, 26 

(n) 

although aesthetic is not very useful since the hard sphere correlations G HS a are complicated functionals of the 
field ip. However the case n = 1 is of some interest. In that case eq. (3.12) says that 

Pfi [K}](1) = (Phs, [{v a + i(t>*}} (l)) H , (3.13) 
It follows readily from the expression (3.4) of the charge correlation function that we also have 

G<?> [K>] (1, ■ • • , n) = (Gg c [{V a + (1, . . . , n)) u , (3.14) 

where 

c (l, . . . , n) = « Q1 . . . q an r ai (1, 1') . . . r Q „ (n, n') ai ___ Qn (l', . . . , n') . (3.15) 
Specializing eq. (3.14) for n = 1 we note that 

PC [K}](1) - (^S, C + i&J] (1)>„ , (3-16) 

where 

PffS, C(l) - ga7V,(l, l'Ws, a(l') • (3-17) 
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2. Relations between G c and G v 



, , , , ,, ^C") nn j 

(n) 



It follows readily from the expression (3.7) of G c and from the KSSHE representation (2.28) of the grand partition 
function that 

^n/2 G (n) _ _ _ = s _x [{ ^ }] ^ | ^ ^ ^_ 1 (^-l^ 

X 6(1) . . . 0(n) S ffS [{^a + #a}] • (3.18) 

At this point we make the remark that 

S E H s [{va + i(j>a}} _ S i(t> a {V) 5 E HS [{V a + i<j> a }] 
6<p(l) ~ 8<p(l) 8v a (V) 

= 6(1) E HS [{v a +i<p a }} . (3.19) 

The relation (3.19) enables us to replace the operators Q(i) which occur the right hand side (RHS) of eq. (3.18) by 
functional derivatives with respect to the field <p. Then, performing n functional integrations by parts 23 ' 26 yields 

^c'c ^-r{ "t^M^ }) } H 

The relation (3.20) can be used to obtain an explicit representation of G^ in terms of the field correlations as long 
as n is not too large. Let us first consider the case n = 1 in which Eq. (3.20) takes the simple form 

Ax < ip(l) > H = -4m0 2 p c {l) . (3.21) 

Once again (cf eq. (2.22)) we obtain the Poisson equation, the solution of which is of course 

<^{l)>n=iP 1/2 V{l), (3.22) 

where V(l) is the GC average of the configurational electric potential, i.e. V(l) =< V(l) >gc- 
In the case n = 2 eq. (3.20) says that 

/?4 2) (1,2) = =^5(1,2) - ^A^Gf (1,2) , (3.23) 

or, by reverting the equation 

GP)(1, 2) = v c (l, 2) - f3v c (l, l')G%\l', 2')v c (2', 2) . (3.24) 

Eqs. (3.23) and (3.24) extend to electrolyte solutions relations obtained recently for neutral fluids. 26 Equations of this 
type were also derived by Ciach and Stell in the framework of a heuristic field theory of the RPM. 29 

By combining eqs. (3.21) and (3.23) one can show easily that the truncated two-body charge correlation function 
satisfies to a similar relation, i.e. 

/?4 2) T (l, 2) = =^5(1, 2) - Al A 2 G< 2 > r (l, 2) . (3.25) 

In the case n = 3 eq. (3.20) yields an awkward expression for Gq\ However the truncated 3-body charge correlation 
function takes the simple form 

if3 3 / 2 G^ T (l, 2, 3) = ^ Al A 2 A 3 G(f) T (l, 2, 3) , (3.26) 

which can be obtained by brute force calculation. The above result suggests that there are simple relations between 
G^ T and G^ T for values of n > 3. Indeed, let us apply the operator 6(2) . . . 6(3) (n > 3) to both sides of eq. 

(3.21). Then, making use of the hierarchy equations satisfied by the G^ T and the G^ T (cf eqs. (Al) and (A9) of 
appendix A) one gets immediately the aesthetic generic formula 

inf3 n/2 G (n) r (1) _ ^ n) = (_^ Al . . . A „ G W T (l, , . . . , n) (Vn > 3) . (3.27) 
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3. Correlations of the electric potential 



It is obvious that 

(V(l) . . . V(n)) cc = v c (l, 1') • • • v c (n, n') G< n) (l', . . . , ri) 



V(l)...V(n)Y =v c (l,l')...v c (n,n')G<g )T (l',...,n'). (3.28) 



T 

= v c (l,l')...v c (n,n')G c 
Combining the above relations with those obtained in sec. (Ill B 2) one gets 



*/? 1/2 (V(V) GC = Ml)>« , (3.29a) 

p(v{l)V{2)) T =v c (l,2)-G^ T (l,2), (3.29b) 

z"/3"/ 2 (y(l)...i>(n)) T =GW T (l,,..,n) (Vn > 3) . (3.29c) 
What can be learned from the above relations is the subject of next section. 

IV. STILLINGER-LOVETT SUM RULES 

A salient property of 3D ionic liquids is the screening effect. To paraphrase Ph. Martin, this type of fluid " in 
thermal equilibrium does not tolerate any charge inhomogeneity over more than a few intermolccular distances". 30 
Even at the liquid- vapor critical point where the correlation length associated with the fluctuations of density diverges 
it is believed, and has been checked by means of numerical simulations in the case of the RPM, 31 that the correlation 
length associated with the fluctuations of charge remains finite. From the existence of screening it is possible to 
deduce sum rules for the charge correlation functions for both homogeneous and inhomogeneous systems, the so- 
called Stillinger-Lovett (SL) sum rules. ' ' As pointed out by B. Jancovici, these rules may be rederived under the 
sole assumption that the system behaves macroscopically as a conductor in the sense that " the laws of macroscopic 
electrostatics are assumed to be obeyed for length scales large compared to the microscopic characteristic lengths of 
the model." 33 We examine below how the SL rules can be deduced from simple hypothesis on the behavior of the 
KSSHE field correlation functions. 

Let us first consider a homogeneous system. In this case < ip >u is a constant as well as are the densities p a 
and the smeared density of charge pc- It follows then from eq. (3.21) that the smeared density pc = 0. This is 
nothing but the usual electroneutrality condition since pc = p a quT a (0) can be identified with the usual local charge 
density for r Q (0) = 1 (property (2.1b) of the smearing function r a ). Therefore, in the framework of our formalism, 
the electroneutrality condition 

PocQcc = (4.1) 

is automatically satisfied for an arbitrary set of chemical potentials {va}, a well known property of ionic fluids. 11 Note 
that eq. (4.1) implies that there are only M — 1 independent chemical potentials. 

The correlations of the electric potential has been studied by various approaches and asymptotically one has 

P(V(1)V(2)) T GC =« c (l,2), (4.2a) 

V(l)...V(n)) Gc =0. (4.2b) 

It must be stressed that these expressions are valid for relative distances large compared to the microscopic charac- 
teristic lengths of the system and if the correlations decay fast enough or, equivalently, if the system behaves as a 
macroscopic conductor. 11 ' 30 ' 33 ' 34 The comparison of these asymptotic behaviors with the exact relations (3.29) derived 
in sec. (Ill B 3) entails that the truncated KSSHE field correlation functions G^ T are short range functions; stated 
otherwise, the KSSHE field is a non-critical field. Conversely, this property being taken as given, we show now that 
one can infer the SL rules. 

We consider now a non-homogeneous system and we take the Laplacian of eq. (3.29b). We get 
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P (p c {l)V{2)) T Gc = 5(1, 2) + i- AxGW t (1, 2) . (4.3) 

Let us integrate eq. (4.3) over f\. With the hypothesis that G v (1, 2) is short range the integration of the Laplacian 
gives zero by an application of Green's theorem. Therefore 

pj d{l) (pc(\)V{2)f Gc = 1 (4.4) 

which is the Carnie-Chan sum rule. 19 In the case of a homogeneous system, the Carnie-Chan sum rule is equivalent 
to the SL sum rule. 33 Let us retrieve these sum rules in our framework. We start with eq. (3.25) for a homogeneous 
system. In this case T (l, 2) = G^ T (r = r\ — r 2 )) and we have 



-AirpG^ T (r) = AS{f) - i-AAG^ T (r) . (4.5) 
With the hypothesis that G^p T (l, 2) is a short range function the integration over r gives zero, i.e. 

J d 3 rG { c )T (r) = (4.6) 
which is the first SL rule. 18 Similarly, after integration by parts 

47T/3 J d 3 rr 2 G { V T (r) = -6 + |- | rf 3 f AG™ * (f) . (4.7) 

If G$ T (r) is a short range function then the integral in the RHS vanishes and we are left with the second SL sum 
rule 

WJ dWG W T (r) = -l. (4.8) 

The SL rules are more conveniently written in terms of the pair distributions h a ^ introduced at eq. (3.2) 

p 7 q 7 h an (k = 0) = — q a (4.9a) 
^ q a q~tP a P~t J d 3 rr 2 h ar( {r) = -1 (4.9b) 

Although eqs. (4.9b) and (4.8) are equivalent, eq. (4.9a) which says that the cloud of charges which surrounds a given 
ion of species a has a total charge — q a implies eq. (4.6) but is more precise. It can nevertheless be derived directly 
with methods similar to those used in this section. The proof is given in appendix B. 

We want to precise that the results derived in this section are not valid for 2D systems which can undergo a 
Kosterlitz-Thouless transition. 35 In the low-temperature KT phase of a 2D PM, the SL rules are violated and the 
sine-Gordon field should exhibit long range correlations with an algebraic decay. In this case, ip is a critical field and 
< ip > is related to the order parameter of the KT transition. 36 

V. THE MEAN FIELD THEORY 

A. Mean field equations 

We define the MF level or saddle point approximation of the theory by the equation 

£mf[{^}] = exp(-W(^)) , (5.1) 
where, at tp = Tp, the action H is stationary. It follows from the expression (2.29) of H that the stationarity condition 

SH _ = (5.2) 



8ip(r) 
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may be rewritten as 

Ai?(l) = -Ainp 1 ' 2 p HS . c [{v a + (1) , (5.3a) 
a (l) = g a T„(l,l')?5(l')- (5-3b) 

Two comments are in place here. Firstly, it can be pointed out that eq. (5.3a) is very similar to the well known 
Poisson-Boltzmann equation. 37 Secondly, it is instructive to compare eq. (5.3a) with the following exact relation, 
easily deduced from eqs. (3.16) and (3.21) 

Ax (<p(l)) n = -Amp 1 / 2 ( PHS . c [{V a + i<f> a }} (l))n ■ (5.4) 

One notes that eq. (5.3a) could have been guessed from the exact eq. (5.4) with the usual assumption of an absence 
of field correlations at the MF level. 



B. The mean field grand canonical free energy 

The MF grand potential is easily obtained by substituting Tp in eq. (5.1) with the result 
log 5 M f[{^}] = logZ HS [{V a + #«}] ~\{ T P K 1 ]^) 

= logZ HS {{v a + i~4> a }] + - (phs, c [{v a + i4> a }] \v c \p H s, c [{><* + i~4> a }]) (5-5) 

The MF number density of species 7 is obtained by taking the functional derivative of log SMF^a}] with respect to 
the local chemical potential v 1 . On the one hand we have 

j—j-j = PHS, 1 [{^ a +l<P a }\(l) + pHS,0[{^ a +t<P a }\( 1 ) g ^ 

= PHS, 7 [{^ + i0 a }](l) + if' 2 PHS, c[{V a + i^ a }](l')Mj| , (5.6) 

and, in the other hand 

Substracting eqs. (5.6) and (5.7) and making use of the stationarity condition (5.3a) gives us 

Pmf, 7 (1) = Phs, ~t[{v a + i~4> a }](l) , (5.8) 
from which we infer the expression of the MF charge density : 

PMF, c(l) = PHS, c[{Va +# C J]( 1 ) ■ (5-9) 

In the case of a homogeneous system the MF KSSHE field Tp is uniform and <j> a reduces to <p a = ^^qaTp since 
t q (0) = 1. Moreover, in this case, the MF equation (5.3a) takes the form of the charge neutrality condition 

Pmf, c = Qa Pmf, a = , (5.10) 

where pmf, a = Phs, a[{v a + i^^qcJp}]- Therefore Tp is an imaginary number which adjusts itself in such a way that 
the charge neutrality condition (5.10) is satisfied. For instance for a RPM with v\ = v<i one has Tp — 0, while for a binary 
SPM (M = 2) one finds, with the same hypothesis on the chemical potentials, that Tp = — i^f]3\og{\q2 / / (q\ — 92)- 
The MF grand canonical free energy Amf is defined as the Legendre transform of logEMF[{v a }} with respect to 
the local chemical potentials -jV Q }. 24 ' 27 ' 38 Therefore one has 

(3A M F =< PMF,aWa > - l°g S MF [{v a }] ■ (5.11) 

When expressed in terms of its natural variables {pmf, a}, the functional (3Amf[{pmf, a}} reads as 
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(3A M f[{pmk a }] = (3Ahs[{pmf, «}]- < Pmf, aWa > + ^(pmf, c\v c \pMF, c) , (5-12) 

as a short calculation will show. The functional Af/s[{/° q }] which appears in the RHS of eq. (5.12) is the exact GC 
free energy of an inhomogeneous hard spheres mixture. As it is well known it is a convex functional of the densities 
{p a }- 38,39 On the other hand the quadratic form < pc\v c \pc > is definite positive, therefore (strictly) convex. The 
last contibrution to Amf, i.e. the linear term < p a \v a >, will not change our conclusion that *4m_f[{/0q}] is a convex 
functional of the densities {p a }- Thence one has, for all < A < 1 

Amf [{Ap a + (1- \)p' a } < A Amf [{p a }} + (1 - A) Amf {p' a } ■ (5.13) 

Since *4m_f[{p q }] is convex then its Legendre transform exists and is also convex. Of course, the Legendre transform 
being involutive, it must be identified with logEMF[{v a }}- Therefore one has 38 

(3Amf[{Pc}] = sup {{p a \v a ) - logS M _F[i> Q }]) , (5.14a) 
{"«} 

logS Mf [{v a }] = sup {{p a \v a ) - (3 Amf [{Pa}]) • (5.14b) 

Some comments are in order. Firstly, assuming that a mixture of hard spheres cannot undergo a liquid- liquid transition 
(liquid-vapor or dimixion) 1 one can safely conclude that Affs[{/°a}] is a strictly convex functional in the fluid region. 
As the sum of two strictly convex functionals (cf eq. (5.12) ) the MF free energy functional Amf[{Pc«}] is also strictly 
convex in the fluid region as is the Poisson-Boltzmann functional. 40 In eq. (5.13) the inequality can thus be replaced 
by a strict inequality for < A < 1. An important consequence is that the solution the the MF equation (5.2) is 
unique for a given set of chemical potential {is a } which rules out the possibility of a fluid- fluid transition of the PM 
at the MF level. Secondly, the charge neutrality condition (5.10) implies that for a homogeneous system the Helmoltz 
free energy per unit volume takes the (too much) simple form 

(3f MF ({ Pa }) = = l3f HS ({ Pa }) Pa u s a . (5.15) 

where the charge neutrality condition p a q a = has been imposed (otherwise PJmf diverges to +oo). 



C. Mean field correlation functions 

1. General case 

It is the place here to recall that a necessary and sufficient condition for the convexity of logS M_F[{f Q }] and 
-4m_f[{Pq}] is that their second order functional derivatives are positive operators, i. e . 24 ' 38 - 39 

SV ° 1 ^(1)^(2) ' S V > ° ' 

spa 1 *, a (i)W) 1 S V > ' (5 - 16) 

We stress that the inequalities (5.16) are strict in the fluid phase because the two functionals logEMF[{v a }} and 
*4m,f[{Ai}] are both strictly convex. The second order derivatives of logSMF with respect to the local chemical 
potentials are the density correlation functions at the MF level: 

Note that, since we consider only two-point functions in this section, we have further economized the notation by 
dropping the subscript (2) in the definition of G MF a/3 . The second order derivatives of (3 Amf with respect to the 
densities are related to the the two-points direct correlation functions. 24,28,38 Define 

^•"'■^- .Si ,5 - i8) 
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In the terminology of statistical field theory Cmf, ap is minus the two-point proper vertex; it is related to the usual 
direct correlation function c a/ g of the theory of liquids by the relation 38 

Cmf, a /s(l, 2) = c MF . af3 (l, 2) - — U<W 5(1, 2) . (5.19) 

Pal, J- J 

The strict convexity of the functionals logSM_F[{^a}] and /3AMF[{p a }] guarantees the existence and unicity of the 
functions G MF ^ Q(3 and Cmf, a/3- Moreover the operators G^ F) af3 and —Cmf, a/3 are both strictly positive and 
-Cmf, a/3 is the inverse of G T MF ^ a0 , i.e. 

-Cmf, a/3(l, 2) G^f, ^ 7 (2, 3) = £> a>7 <5(1, 3) , (5.20) 

which, when reexpressed in terms of the functions h a p and c Qi g, takes the familiar form of the Orstein-Zernike equation 1 

h,MF, 0/3(1, 2) = cmf, q/3(1, 2) + /imf, o 7 (1, 3)p 7 (3)cmf, 7/3(3, 2) . (5.21) 

Cmf, 0/3(1, 2) is obtained readily from the expression (5.12) of the MF free energy functional (3Amf [{pa}] with the 
simple result 

Cmf, 0/3(1, 2) = C HS , 0/3(1, 2) - f3q a qpw a0 (l, 2) , (5.22) 

where Crs. 0/3(1,2) is minus the exact two-point proper vertex of the reference HS fluid at the mean field densities 
{pMF,a}- Eq. (5.22) implies that for a homogeneous system and for r\i > (a a +073) the direct MF correlation 
functions reads as 

cmf, 0/3(7-12) = c HS , 0/3(7-12) - Pq a qp/ri2 , (5.23) 
which is nothing but the RPA closure of the theory of liquids. 1 ' 21 

iMF ~ 



The calculation of G^ F ' T (l, 2) is more involved. Our starting point will be the equation 



SpMF, o[{> 7 }](1) 



GW, a/9 (l,2)- ^ (2) 



<M 2 ) 

Clearly pns, a [{^7 + *0 7 }] depends upon z/g directly but also through the smeared mean fields (j> 7 ; therefore one has 



Cmf, 0/3(1^ 2) — 



SpHS, a[i> 7 + &<M](1) 



<M 2 ) 



fyffs, o[{^ 7 +i(/i 7 }](l) 



my 



<M2) 



G^ SiQ/3 (l,2)+i/3 1 / 2 G^ Fia7 (l,4) g7 r 7 (4,3)^^ . (5.25) 



Now, remarking that the stationarity condition (5.3a) implies that 

= ^/ 3 1 /2 ^ , c (3, 5)« 4 t 4 (5, 6)G^ F] Jj8 (6, 2) , (5.26) 



one obtains 

Gm_f, 0/3(1, 2) = G^ s a/3 (l, 2) 

- Q7 (l, 4)<z 7 r 7 (4, 3)w c (3, 5)^(5, 6)G^ Fi Jj8 (6, 2) , (5.27) 

which can be rewritten under a matricial form as 

G^ F (1, 2) = G T HS (l, 2) -13 Gl s (l, 3) v(3, 4) G^ F (4, 2) , (5.28) 

where G^ /F ( H5 )(1, 2) denotes the matrix of elements Gj^p^^ Q/3 (l, 2) and w(l, 2) that of elements v a p(l, 2) (cf eq. 
(2.4)) The formal solution of eq. (5.28) is then 
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G T MF (1,2) = (l + 0G T HS *v) 1 * G T HS , (5.29) 

where l a/3 (l,2) = S a p5(l,2) is the unit operator and the star "*" denotes a convolution. One easily checks that 
GJ^p a g and Cmf, a/3 as given by eqs. (5.29) and (5.22) respectively do satisfy the OZ equation (5.20). 

In the case of a homogeneous system cq. (5.29) can be simplified considerably. Let us work in Fourier space where 
we have 

G T MF (k)= (U-P^y 1 .G T HS (k), (5.30) 

where U_ aj3 — 8 a p is the unit matrix of rank M x M and P(k) denotes the matrix of elements 

Pup = -Pvc{k)&fj St al (k)q 1 q f3 T^(k)T f 3{k) (no sum over (3) . (5.31) 

P(k) has the remarkable property that 

E(k) = -{iv c {k)G T HS , c (k)P(k) , (5.32) 

which leads us to search the inverse of U_ — P_{k) under the form (U_ — P(fc)) -1 = U + Jj.(k)P_(k). The identity 
(U + Ji(k)P_(k)).(U_ — P_{k)) = 1 when combined with the property (5.32) implies that 

H{k) = L- . (5.33) 

l + f3v c (k)G T HS ^ c (k) 



Therefore 



(u-P(k)) =U+ =U- , (5.34) 

from which a simple expression for G]^ F a p(k) is easily obtained 



r T <k)-r T (k) Pvc{k)T a {k)T p (k) 

l + 0v c (k)G% SiC (k) 



T a {k) = G 1 HS ,^{k)q 1 T 1 {k). (5.35b) 

The expressions of the Fourier transforms of the usual pair correlations Iimf, a[) can then be deduced from eqs. (5.35). 
One finds 

hMF, a p(k) = h H s, a 0{k) - — . 5.36) 

Pa pp 1 + f3v c {k)G T HSC (k) 

An important comment is in place here. As well known, the fact that the direct correlation functions (5.23) behave 
as the Coulomb potential at large distances is sufficient to ensure that the SL rules are satisfied. 1 ' 41 Therefore, the 
Hmf. a/3(j), the Fourier transforms of which are given above, automatically satisfy to the SL sum rules (4.9). 

The two-point MF charge correlation is obtained by taking the convolution of the two members of eq. (5.27) with 
q a Q0T a (l, l')rg(2, 2') which gives the equation 

G T MF , c (l, 2) = G T HS , c (l, 2) - pGF HS< c (l, 1') v c (l', 2') G T MF c (2', 2) , (5.37) 
the formal solution of which is 

G T MF , c = (1 + 0G T HS , c * «c) ~* * G T HS , c . (5.38) 
Therefore, for a homogeneous system, the Fourier transform of Gj 4F c (r) has the simple expression 

G T MF c (k) = fc G hs c(k) 

MF ' C[) fc 2 + 4 7 r/3G^ c (fc) 
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Finally the truncated pair correlation function of the KSSHE field is obtained by using eq. (3.25) which gives us 

G T MF , 2) = v c {l, 2) - /3« c (l, 1') G T MK c (l', 2') v c (2', 2) 

= (1 + /JG£ Si c * Uc ) _1 (1, 1') v c {V, 2) . (5.40) 
For an homogeneous fluid, one thus has in Fourier space 

4-7T 

G T MF , v {k) = ~ T — • (5.41) 

2. Application to the special primitive model 

In the SPM all the spheres have the same diameter a a — a and we also assume that all the smearing distributions 
T a are identical with the same r. With these assumptions it is shown in app. D that the functions G^ S T C are nearly 
decoupled of the density correlations G^g T of the HS reference fluid. For instance one has for n = 2 

Gls,c(k)=PMF, a q 2 a ^(k) , (5.42) 
which entails considerable simplifications. Indeed the MF pair correlation functions take the simple form 

hMF, af }{k) = h HS {k) - fc2 + ^ , (5.43) 

where k mf = 4irf3pMF, aQa ^ s the squared Debye number and hHs( r ) the truncated pair correlation function of the 
HS fluid at the density Pmf — ^,PMF,a- Eq. (5.43) when combined with the electroneutrality condition implies 
that the correlations of the (total) density are equal to those of the HS reference fluid which rules out a liquid-vapor 
transition of the SPM at the MF level. Moreover the charge and KSSHE field correlations of the SPM appear to be 
completely decoupled from the density fluctuations since one finds 



47T/3 k 2 +? 2 (k)K 2 MF 



4"7T 



Note that for point charge distributions (i.e. f (k) = 1) the above functions reduce to the Yukawa potential, more 
precisely 

g mf, v ( r ) = exp(-K M F r)/r , 

and 

4tt/3Gmf, c( r ) = K iiF^ 3) (r) ~ K 4 MF exp(-K M F r)j^r . 



D. The mean field free energy as an exact lower bound of the free energy 

It can be shown, under quite general conditions, that the logarithm of the grand-partition function logS[{/x Q }] is 
a convex functional of the local chemical potentials even before the passage to the thermodynamic limit. Similarly, 
the exact Kohn-Sham free energy /3„4[{p Q }] is a convex functional of the local densities. Moreover log and 
-4[{/0 Q }] constitute a pair of Legendre transforms which can be expressed as 24 > 38 > 39 

l3A[{pa\] = sup {(p a \v a ) - logS[{^ Q }]) , (5.45a) 
logS[K}] = sup ({p a \v a ) -0A[{p a }]) . (5.45b) 

{Po} 
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Recall that the Young inequalities which follow directly from eqs (5.45a), (5.45b) say that 

pA[{p a }] + logS[K}] > (p a \is a ) (VK>, V{ Pa }) . (5.46) 

We have seen in sec. (VB) that the MF functionals log Smf [{Ma}] an d Amf[{p q }] are also two convex functionals 
linked by a Legendre transform. We establish below some rigorous inequalities between the exact and MF free energies 
and grand potentials which extend to ionic fluids results recently obtained for neutral fluids. 26 
First, we notice that the fundamental eq. (2.28) rewritten as 



S[K}] - (exp (log S HS [{v a + i4 +*4>a}]) ) Vc (VK» , (5-47) 
can be inverted in order to give 

H/fs[K}] = (cxp (log ~[K - v s a + K}])) Vc (V{i/ a » • (5-48) 

Then, we apply Young inequality (5.46) to the RHS of eq. (5.48) which yields 

E H s[{v a }} > exp {-/3A[{p a }] + (p a \v a - v s a )) x 

x (cxp(0»))„ c (V{i/«}, V{p Q }), (5.49) 

where it follows from the definition (2.17) of the smeared field 4> a that 

i(l)=/3 1/2 ? a r a (l,l'K(l') • (5.50) 
The field average in the RHS of eq. (5.49) is computed by making use of Wick's theorem which yields 

(exp (0M))„ c = exp (pc\v c \pc)j ■ (5.51) 

By inserting the above result into eq. (5.49) we obtain 



logZ HS [{v a }] > -0A[{p a }] + (Pa\"a - > + f (pMpc) (VK>, V{ Pq }) , (5.52) 



which implies 



f3A[{p a }} > sup \{p a \v a ) - logins [{Va}]] 
{"«} 

+ {(Pc\v c \pc)-( Pa \v S a ) (V{p a })- (5.53) 

The sup in the RHS of eq. (5.53) is the Legendre transform of the exact grand potential logSfl-g[{i/ Q }], i.e. the free 
energy /3A HS [{p a }]. It follows then from the expression (5.12) of /3Amf that eq. (5.53) may be recast under the form 

0A[{pa}]>0AMF[{Pa}] (V{Pa» ■ (5-54) 

Therefore the MF free energy /3*4MF[{Pa}] is a rigorous lower bound of the exact GC free energy /?^1[{/3 Q }]. Moreover 
since /3Amf and logSps are linked by a Legendre transform (see eqs. (5.14a) and (5.14b)) we also have 

log 2 M F [{>«}] = sup [(p a \v a ) - l3A M F[{p a }]] (V{^ q }) 

> sup [(p a \v a ) - PA [{p a }}} (V{i/ a » , (5.55) 

which can be, with help of eq. (5.45b), rewritten as 

log£ M F[K}]>logS[K}] (VK». (5.56) 

In other words logS#s[{f Q }] is a rigorous upper bound of the grand potential logS[{f Q }]. Recall that, for a homo- 
geneous system in the point-like limit, the Debyc-Hiickcl approximation for the pressure constitutes a rigorous lower 
bound for the pressure. This result can also be proved with the help of the Sine-Gordon transform. 42 
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Let us specialize now to the homogeneous case. We denote by f3f({p a }) the exact free energy per unit volume of 
the system. It is bounded from below by /3fMF{{Pa}) and we have therefore 

PfdPa}) > P.fHs{{p a }) - Pa v S a , (5.57) 

where the electroneutrality condition pc = has been imposed. It follows from cq. (5.57) that the "best" or 
"optimized" MF free energy is obtained by minimizing the functional K [{r a }] = p a v^ — j3p a q 2 t w aa (0)/2 with respect 
to the variations of the smearing functions {r Q (r)}. Since these distributions are normalized to unity we introduce 
M Legendre parameters A Q and define the functional 



Q, 



t [r a ] = f d 3 r Ta (r) . (5.58) 

Jr<a a 



Now, we minimize i4T[{r Q }] for a given set of X a . The constraints Q a = 1 (Va ) will utimately serve to determine the 
values of the Legendre parameters \ a . The conditions for an extremum of K are 



SK[{T a }} 



Since w aa (ri 2 ) = w aa (l,2) = r Q (l, l')r Q (2, 2')v c (l' , 2') we have 

6w aa {0) 



(V 7 G (1,...,M), \\r\\ <a 7 ) . (5.59) 



28 ai V 7 {r) , (5.60) 



where V 1 is the electrostatic potential created by the charge distribution r 7 (i.e. Ky(r 12 ) = V 7 (l, 2) — r 7 (l, 2')v c (2', 2)). 
The stationarity condition (5.59) takes thus the form 

y 7 (f) = A 7 (V 7 G (1, . . . , M), ||r1 < a 7 ) . (5.61) 

The potential V 7 created by the distribution t 7 must therefore be constant inside the sphere of radius a a . From 
elementary electrostatics we conclude that r 7 (r) must be an uniform surface distribution of charge of radius a a . In 
order to satisfy the constraint Q~ f = 1 one must have 

r 7 (f) = J(|M|-a 7 )-^. (5.62) 

Note that the solution (5.62) is indeed a minimum of K since its second order functional derivative with respect to 
r a is the Coulomb potential v c (l,2) which is a positive operator. The simple form of the Fourier transform of r 7 (r) 
which is found to be 

r 7 (*) = ^ (5.63) 

allows us to compute explicitly the pair potentials w a p(r) as the inverse Fourier transforms of the functions w a p(k) = 
4TTT a (k)Tp(k) I k 2 with the result 

r > G a + a i3 => w afj (r) = - (5.64a) 

- ^ ^- , - ^ ( \ Va+Vp-r 1 trcAU\ 

<?a - 073 < r < a a + (70 w a p{r) = h (5.64b) 

G a U(i 2r 4r a a crp 

2 

< r < a a — a (3 => w a p{r) = — , (5.64c) 

o a 

where it was assumed that a a > a (3. As a subproduct of these equations we obtain the self-energies — flq^/aa and 
thus the expression of the optimized MF free energy 

PfMF({ Pa }) = [3f HS ({p a }) - ^ Pa ■ (5.65) 

On 
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Note that the corresponding excess internal energy per unit volume /3umf — <9(/3/mf)/<9/3 = —(3p a q\j o a coincides 
with the Onsager lower bound. 43 We note that the MF approximation is thermodynamically inconsistent in the sense 
that the MF energy umf is obviously not equal to that which can be obtained by the integration of gMF, af3(f)w a /s(r). 

Reporting the expressions (5.64) of the pair potentials w a p(r) in the equations (5.18) one obtains the MF direct 
correlation functions. We shall discuss only the case of the SPM where all the radii a a are equal to the same a. In 
this case one finds 

r > a => c M f, ap(r) = c HS {r) - /3q a q/3/r (5.66a) 

2(7 — r 

< r < a => c M F, a p(r) = c H s(r) - (3q a qp — , (5.66b) 

where ensir) is the exact direct correlation function of the HS fluid at the density p = J2 a p a . The expressions (5.66) 
of the Cap are very similar to that obtained in the framework of the MSA approximation i.e. 41,44 

r > a cms A, ap(r) = c^W - Pq a qfi/r (5.67a) 
< r < a cms A, ap(r) = c£s( r ) - f3q a q p — , (5.67b) 

where c|^(r) denotes the direct correlation function of the HS fluid in the Percus-Yevick approximation which is an 
excellent approximation of the exact chs{t) especially at low or moderate densities. 1 The main difference between 
eqs. (5.66) and. (5.67) is the occurrence of a parameter B in the MSA solution. Clearly the MF and MSA solutions 
coincide for B — 1. Since we have 41 ' 44 

K 2 er 2 

B = (na f + kg - no{\ + 2kct) 1 /2 ' (5 ' 68) 

where k is the Debye number, it happens only in the limit k — > oo. For a finite k, B(k) is positive and comprised 
between and 1. Obviously the electrostatic contribution to c^ A (r) in the core (i.e. the second term in the RHS 
of eq. (5.67b)) may be seen as the interaction energy w a p{r) of two surface distributions of charge of equal radii 
a 1 = a IB. However, since B < 1 => a' > a, this interaction energy is not equal to 1/r outside the core. Consequently 
the MSA solution cannot be interpreted as a MF KSSHE theory except in the limiting case B = 1. 



VI. THE GAUSSIAN APPROXIMATION 



A. The general case 

Let us define the Gaussian approximation of the KSSHE theory in the following way. 20 We write <p = Tp + Sip where 
Tp is the mean field solution and Sip a real scalar field and we expand functionally the action H[ip] (cf eq. (2.29) ) up to 
second order in Sip around the MF solution. In this way the exact action H[ip] is replaced by an approximate action 
Hg[<p] given by 

H G [ip]=nm + ^{Sip\A- 1 \Sip) , (6.1) 

where the terms linear in Sip are absent as a consequence of the stationarity condition (5.2), and the inverse of the 
propagator A is given by: 

A- 1 (l,2) = V - 1 (l,2)+^ SiC (l ) 2). (6.2) 
A comparison of eq. (6.2) with eq. (5.40) yields the expected result 20 

A(l,2) = G^ FjV (l,2). (6.3) 
In the Gaussian approximation, the grand partition function is therefore given by 

■=G [{V a }\ = CX P {-HYP]) -j— 7 - x n • ( 6 - 4 ) 

J Vip exp {-±(p\v c \ip)) 
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From now we specialize to the case of a homogeneous system for which an explicit expression of Eq can be obtained, 



i.e. 



3 G [K}] = exp (-«[?]) (6.5) 

where the normalization constants A/a and Af Vc are defined according to eq. (2.13). A well known property of functional 
integrals of Gaussian functional gives us 20 > 22 . 23 



= -^|d v fclog(l + ^G^ c (fc)) . (6.6) 

Taking the logarithm of cq. (6.5) and taking into account the relation (6.6) we thus obtain the pressure of the ionic 
solution at the Gaussian level: 

d v Hog (l + ^q a qf,T a (k)T P {k)G% St Q/3 (fc)) , (6.7) 

where we have expanded the HS charge correlation function in oder to make explicit the dependence of the result upon 
the smearing functions r a (k). Note that the HS truncated pair correlation functions Gjj S a p(k) must be evaluated 
at the chemical potentials {v a + ifi 1 / 2 ^} . As a short examination reveals, the integral in the RHS of eq. (6.7) is 
convergent at large k for any reasonable smearing function r Q (i.e. surface, volume distributions etc) but diverges for 
point charge distributions (r a (fc) = 1). 

Some words on the field correlation functions at the Gaussian level. Following Ma 20 we define 

which entails the relations 

(6<p)n a = , (6.9a) 

C%] J(l,2) = A(l,2) ee G^/ V (l,2) , (6.9b) 

4" ) v T (lv,i)=0forn>3. (6.9c) 

In particular, eq. (6.9a) says that the number densities in the Gaussian approximation coincide with their MF values. 
The Legendre transform of (3Pq is therefore easily obtained and reads as 

PfG [{ P a}} = PfHS [{Pa}} - Pa^ + \ J d" 'k log (l + c (fc)) , (6.10) 

where the truncated charge correlation function G^g c of the reference HS fluid must be evaluated at the densities 
{p a }- Replacing the self energy by its expression (2.6), [3Jq can be recast under the form 

PfaiiPa}} = PfHs[{pa}] + -^PaPpQaQp / d 3 f ll HS , a p{r)w a p(r) 

+ \Jd-k {log (l + ^-Gls, c (fc)) - ^-G T HSi c (fc)} , (6.11) 
which will be discussed in next section. 



A(fc) 
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B. Comparison with the RPA 



As in the case of neutral fluids, 26 the expression (6.11) of the free energy at the Gaussian level coincides with that 
obtained in the framework of the random phase approximation (RPA) of the theory of liquids. The proof of this 
equivalence is however more tricky in the present case. The RPA theory can be summarized as follows. 1,21,37 The 
connected two-body correlation functions are given by 

GW a p0; 2) = G T HSj a0 (l, 2) + PaPa C a0 (l, 2) , (6.12) 

where the "renormalizcd" potentials C Q/ a(l,2) are given in Fourier space by 

-1 ~T 



p a p a C aP {k) = P(k) . (U - P(fc)) . G T HS (k) 



(6.13) 

a/3 



where the matrix P(k) the elements of which are given by 

Pafi = -^s, a p{k)v aP (k) (6.14) 

turns out to coincide with that we defined at eq. (5.31). Note that, in the diagrammatic formulation of Chandler 
and Andersen, G^ s is the matrix of the so-called "hypervertices" . 21 It is easy to check that the properties (5.32) and 
(5.34) of the matrix P(k) imply that 



P(k) ■ (U-P(k)) 1 .G HS {k) 



a/3 l+(3v c (k)G T HS n{k) 



! P(k).G T HS (k) 



a/3 



pv c {k)T a {k)T {k) 

and therefore, cf eq. (5.35), that 



l + (ft c (k)G T HSt c (k) 



(6.15) 



G«PA,a/j(1.2) = G^ iaj8 (l,2). (6.16) 

The expression found by Chandler and Andersen for the RPA free energy reads as 21,37 

a r 

PfRPA [{Pa}] = PJhs [{Pa}} + / d z rh HS . a fj(r)v afj (r) 

+ ^Jd v k { trP(fe) +logdct (t/-P(fc))} • (6-17) 
We first note that it follows from the definition (5.31) of P(k) that 

trP(fc) = -ffi c (k)(f HS} c(k) . (6.18) 

The determinant in the RHS of eq. (6.17) can be explicitly computed from a classical formula of linear algebra. 45 
Recall that for a n x n matrix a we have 

det (a i3 + xSij ) = x n + S^- 1 + . . . + S n x + x n , (6.19) 

where Sk denotes the sum of the minors of order k, i.e. 

S k = J2 det K^) • ( 6,2 °) 

qi<...<q k 

When applied to the matrix Pik) (with x = 1) the formula (6.19) yields the simple result that 

det (u_ - £(fc)) = 1 - trP(Jfe) . (6.21) 



20 



Indeed, since P a f)(k) is of the form s a (k)tp(k) then all the minors in eq. (6.19) vanish except S\ which identifies with 

the trace of the matrix P Ql a(fc). It follows then from the expression (6.17) and the intermediate results (6.18) and 
(6.21) that the RPA free energy is identical to the Gaussian KSSHE free energy: 

PJrpa [{ Pa }} = (3 fa [{p a }} ■ (6.22) 

It must be stressed that in the RPA the pair potential v a p between two ions is arbitrary in the core (i.e. for 
< r < a a + ap). The simple forms that we have obtained obtained for the pair correlations and the free energy 
(cf eqs. (5.35) and (6.10)) are valid only if v a p may be identified with the interaction of two smeared distributions of 
charge. In that case, the relevant mathematical property which yields important simplifications is that v a p(k) may 
be written as v a p(k) = Pv c (k)q a T a (k)qpTf)(k) for all the pairs (a,/?) , i.e. as the product of two functions s a (k) and 

S/}(k). 

In the optimized RPA (ORPA) the pair potentials in the core are chosen in such a way to ensure that the radial 
pair correlations goRPA, a(j{ r ) vanish for < r < a a +ap. 21 Chandler and Andersen have shown that this condition 
corresponds to an extremum of PJrpa considered as a functional of the core potentials. There is no analytical solution 
to this variational problem as far as the author knows, and the solution must be seeked numerically. It is however 
quite certain that this solution cannot be interpreted as the interaction of two smeared charges as suggested by the 
discussion on the MSA integral equation of sec. (VD) (the MSA is identical with the ORPA if the HS reference fluid is 
described in the framework of the PY theory). Therefore, with smeared electrostatic potentials one can never ensure 
that the go, a ( r ) vanish in the core. 

C. The SPM 

We conclude this section by specializing to the case of the SPM. In this model all the hard spheres have the same 
diameter a and all the smearing distributions r a are equal to the same r. It follows from the expression (5.42) of 

(2) T 

G HS c that the Gaussian free energy of the SPM takes then the very simple form : 

Pfa [{ Pa }] = PIhs [{ Pa }\ + ljd v k jlog (l + ^r 2 (fc)) - ^Jr 2 (fc)| . (6.23) 

Several comments on eq. (6.23) are in order at this point. 

(a) As it is well known, 37 and as a direct calculation will show, eq. (6.23) reduces to the Debye-Hiickel free energy 
PJdh — Pfns — k 3 /127t in the limit of point-like distributions (i.e. when r(fc) = 1 for all k). Recall that in the 
point-like limit, the MF free energy diverges. Therefore the first order term in the loopwise expansion of the free 
energy makes the final result finite. 

(b) The expression (6.23) shows no sign of a hidden "RPA catastrophe" 46 which would be the case if the argument 
of the logarithm happened to be negative for some k. This happy circumstance is a consequence of the regularization 
by "smearing" of the Coulomb potential which was adopted in this paper. Other types of regularization can however 
lead to such a RPA catastrophe. For instance, adopting the Weeks- Chandler- Andersen (WCA) recipe 47 one could 
set the potential constant in the core, i.e. v a p(r) = q a qp/<r for < r < a. Clearly, one has in fact v a f}(r) = 
q a qpw(r) (Vr) where w(r) can be interpreted as the electric potential created by the surface distribution of charge 
ts(t) = S(r — lr)/'rr<T 2 . Thus w(l,2) = rs(l, l')v c (V, 2). As a consequence, one must replace the term r 2 (fc) by 
Ts(k) — sm(ka) / ka (not squared) in eq. (6.23), with the annoying consequence that the argument of the log can 
become negative for some k yielding a "RPA catastrophe". The mathematical origin of this catastrophe is that, 
within the WCA framework, the energy is not a positive definite quadratic form and the KSSHE transform is ill- 
defined. It has been suggested that this instability of the RPA theory could possibly be related to the order-disorder 
transition of the lattice-version of the RPM. 29 ' 48 Note that this transition is not present in the continuous version 
of the model although the RPA catastrophe for the WCA-RPA free energy is still there. In all cases, one cannot be 
comfortable with a theory the ability of which to predict a possible transition depends on the way the regularization 
of the interaction is performed. 

VII. A TWO-LOOP ORDER CALCULATION 

A. The cumulant expansion 

We adopt the same decomposition <p = Tp + Sip of the KSSHE field as that considered in the section on the Gaussian 
approximation but, this time, we do not truncate the action at the second order in Sip and consider rather the full 
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functional Taylor expansion of H[<p] about the saddle point 



H[<p] = H[Tp] + ±(6<p\ AZ 1 !^) + *H[5<p\ , (7.1a) 
AH[8<p] = V -fW (n) [?5](l, ■ • • , n)ty(l) • • ■ M«) , (7.1b) 

n— 3 

where the free propagator defined at eq. (6.2) must be evaluated at the saddle point Tp and the integral kernels 
are given by 

H (n) m(l,...,n) = -(il3 1 ' 2 ) n q ai ...q an T ai (1, l') . . . r a „ (n, n') x 
4I,U [{*« + (l',...,n'), 

^-^rc^^a,...^). (7.2) 

With these notations, the grand partition function can be recast under the form 

S \{u a }\ = exp (-H [$>]) ^ (cxp (-AWN))a_ , (7.3) 

yielding the cumulant expansion 20 ' 49 ' 50 

logS = -n[Tp] + log ^ -(AW) Af +1(AW 2 )^ -^(aw 3 )X f +•■• (7-4) 

The two first terms in the RHS of eq. (7.4) correspond the MF and Gaussian approximations respectively, the other 
terms are Gaussian averages which can be computed with the help of Wick's theorem. Reordering these additional 
terms in powers of some small parameter yields various types of expansions. Low fugacity and high temperature 
expansions of the pressure and the free energy of the RPM where obtained in this way in papers I and II, yielding 
to a rediscovery of old results derived years ago by Mayer graphs expansion technics. 13 ' 14 In both cases the "small 
parameter" has a physical origin and thus each order of the expansion should be independent of the pair potentials 
in the core and therefore of the smearing functions r a . This point was checked carefully in I and II for the RPM by 
considering surface distributions of radius a < a and by showing that the first terms of the fugacity and temperature 
expansions were indeed independent of a. 

Here we consider a loopwise expansion. The small parameter A of this expansion is defined by rewriting 49 ' 50 

Sip -» X^Hip , 

logS -> A logS. (7.5) 

An expansion of eq. (7.4) in powers of A yields the loopwise expansion of logS; at the end of the calculation one set 
A = 1. The "small" parameter A(= fi !!) is not a physical parameter but serves mainly to keep track of different classes 
of Feynman diagrams; therefore each term of the loop expansion can a priori depend on the smearing functions T a . 
One can only hope that the larger the number of terms is retained in the expansion the "smaller" this dependence 
will be. In a sense, the fact that the MF free energy diverges for point charge distributions but remains finite in 
the Gaussian (one-loop) approximation confirms this hope. Henceforth we shall retain for r a the distributions (5.63) 
which lead to the best MF free energy. Finally, we note that it follows from eqs. (2.28), (2.17) and (7.5) that a 
high-temperature expansion of the loop-expansion of log S at order \ k obviously yields the correct high-temperature 
expansion of log S at order j3 k . 

B. The two-loop expansion of log H 

The loop expansion of log S for the general Hamiltonian (7.1) can be found in the literature. 49 ' 50 Taking into account 
the form (7.2) of the kernels Ti^ one finds, at the two-loop order 



22 



A/a— 

logS[{j/ Q }] = logS M F[{t'a}]+Alog-^ L 

A 2 /? 2 (4)T 



G^; c [^](l, 2, 3, 4) A^(l, 2) A^(3, 4) (= D a ) 



G^; c [^(x 1 ,x 2 ,x 3 )G^ T c [^(y 1 ,r 2 ,y 3 ) x 

x A^(X 1 ,X 2 )A^(y 1 ,F 2 )A^(X3,F 3 ) (- A,) 
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x A ¥ (X lt Fi)A^(X 2 , F 2 )A^(X 3 , y 3 ) (= D c ) 

+ G(A 3 ), (7.6) 
where the diagrams corresponding to the two-loop order terms of the RHS of eq. (7.6) have been sketched in fig. (1). 






D b 



D, 



FIG. 1. Two-loop diagrams contributing to log^,. The lines represent the bare propagator A„ 



Some remarks. In general the loop-expansion (7.6) does not converge notably in the critical region as a consequence 
of the singular behavior of the correlation functions. 49,50 Here however, we deal with a non critical field ip the correlation 
functions of which were shown to be short range in section IV; eq. (7.6) is thus expected to be an asymptotic expansion. 
The second remark is that the expression (7.6) is horribly intricate in the general case since, a priori, it involves 
the 2, 3-body and 4-body density correlation functions of the reference HS fluid. However it turns out that eq. (7.6) 
can be considerably simplified for the SPM since, in this case, only the 2-body correlation functions survive, yielding 
tractable expressions. 

Before specializing to the case of the SPM let us introduce some notations. All the thermodynamic quantities of 
interest will be expanded in powers of A. We shall note, for intance for the density of species a 



(7.7) 



with, of course = pmf, a = Phs, a[{v a + iq a p}}- As it is well known, the two-loop order expansion of the free 
energy requires the densities at the one-loop order only. 49 ' 50 We have checked in the general case that the one-loop 
order correction to the densities preserves the electroneutrality; the detailed calculation is made in appendix (C). 

Henceforth we consider only the homogeneous SPM. All the ions have the same diameter a and all the smearing 
functions are equal to r{k) — sin(fc(j) /ka. We show in app. D that the charge neutrality condition at the MF level 
(pQ^ = 0) kills most of the terms involved in the expressions of G^\ T C [p] and that the only surviving terms are 

G^ i T [^](l,2) = pMr(l,l')r(2,l') (= D 2 ) 
G^ s T c [p](l, 2, 3) = p^ql r(l, l')r(2, l')r(3, 1') (= D 3 ) 
G^ T C [^](1,2,3,4) = pWqt r(l,l')r(2,l')r(3,l>(4,l') (= D 4a ) 

+ (pi 0) <) 2 x{ 

r (l,l')r(2,l')/i ff s(l',2')r(2',4)T(2',3) (= D 4b ) + 
t (1, l')r(3, l')h HS (l', 2')r(2', 4)r(2', 2) (= D ic ) + 

r (1, l')r(4, l')h HS (l', 2')r(2', 2)r(2', 3) (= D id )} , (7.8) 



where hns ( r ) = 9hs ( r ) — 1 is the pair distribution of the reference HS fluid at the total MF density //°) = J2 a p, 
Diagrammatic representations of the various terms of the RHS of eq. (7.8) have been sketched in fig. (2). 



(o) 
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D 4 b D 4c D 4d 

FIG. 2. Diagrams representing G^' S T C (n = 2,3,4). Dashed lines represent /i W s(l,2) and solid lines the smearing function 
r(l,2). 



Replacing either algebraically or graphically the kernels G^g T c [if] by their expressions (7.8) in cq. (7.6) yields the 
expression for the pressure (3P = log£/fi at the two-loop order, i.e. P = P(°) + + \ 2 P^ +C(A 3 ) with of course 



and 



0PM [{u a }] = pP MF [{u a }] , 

pp^[{u a }] = - 1 - [ t i i,- ;„ufi • ■ 



(7.9) 



k 2 



(7.10) 



( where 2 = n 2 MF = 4ir/3p a q 2 is the squared Debye number at the 0-loop order) as we already know, and 



+ 



[a(°) ( o) 

2 

/3 3 



d 3 rh HS (r) (=D 2 ) 



'J d 3 rh HS (r)[A?\r)\ 2 } (= £> 3 ) } 



/3 3 r jo) 3 i2 



n 2 
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ft ?o 



A(°)( ) A(°)(0) (=D 4 ) 



A(°)( r) l (=£,5), 



where the "smeared" propagator A^ ' which enters the RHS of eq. (7.11) is given by the convolution 



or, in Fourier space 



A(°)(l,2)=r(l,l')A^(l',2')r(2',2) 
47rf 2 (fc) 



A(°)(fc) 



fc2 + K (0)2 f 2( fc ) • 



(7.11) 



(7.12) 



(7.13) 



We have already met the function A?'(r) (under the name X T (r)) in papers I and II. Recall its high temperature, or 
low behavior 



r < a => A|°)(r) = - - 4 " « (0) + <?(« (0) 2 ) 
a a 

a AW(r) = sinh V 0) g) e*P(-" ( °M x n , , ( .n 2 , 



r > 



( K (°)<t) 2 



(7.14a) 
(7.14b) 



Diagrammatic representations of the various terms of the RHS of cq. (7.11) have been sketched in fig. (3). 
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Dj D 2 




D 3 D 4 D 5 

FIG. 3. Two-loop diagrams contributing to the pressure (cf eq. (7.11)). The smeared propagator Ar°'(l,2) has been repre- 
sented by a solid line plus a slash. As for the diagrams of fig. (2) the dashed lines represent tiHs(l, 2). 



C. The two-loop expansion of the free energy 



As it is well known from statistical field theory, the reducible diagram Db of fig. (1) should disappear after a 
Legendre transform of the pressure yielding a simple expression for the free energy. 49 ' 50 However in statistical field 
theory one is interested with the Legendre transform of f3P with respect to the mean value of the field < ip >t-l\ by 
contrast, in our case, the quantity of interest is rather the Hclmoltz free energy f3f defined as the Legendre Transform 
of the pressure with respect to the M chemical potentials v a (cf eqs. (5.45a), (5.45b)). These two Legendre transforms 
do not coincide a priori and a direct calculation of f3f at the two-loop order is therefore necessary. 

We have first to compute the densities at the one-loop order since the terms X 2 pi" 1 and of higher order do not 
contribute to (3f at order 0(A 2 ) (included) as a consequence of the stationarity condition (this point will emerge in 
subsequent developments). Eqs. (7.9) and (7.10) entail that 



dpP [{!/«}] 



= P (0) - 



AAi o) (0) 9«(°) 2 



0(A 2 ), 



2 4tt dv a 

where we recall that p° a — pHS, a[{pa + i^^qaTp}] is the 0-loop order (or MF) density of species a. We note that 

g K (0) 2 



(7.15) 



dv a 



= 4^/3 q 2 (dp^/dv a ) = 47T/3 q 2 G T MFy aj (k = 0) , 



(7.16) 



where Gjf F ai can be computed from cq. (5.43). After some algebra one finds finally 



(1) _/3M0) 



h HS (0) pM+pgWo 



(7.17) 



(no summation over a in the RHS). We have already pointed out that the clcctroneutrality is satisfied at the MF level 
(i.e. pa\ a = 0) and we note with satisfaction that eq. (7.17) implies that it is also satisfied at the one-loop order 
(i.e. Pa^Qa — 0; see also appendix C for a thorough treatment in the general case). As a subproduct of cq. (7.17) we 
obtain the one-loop order corrections to the total density and Debye number 



n 2 W ee A*p£>& 

= -27r/3 2 AW(0) 



(l + p^h HS (0)) 



h HS (0) + 



(7.18a) 



(7.18b) 
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The free energy at the two-loop order is obtained by reexpressing 

13 f = p a v a - I3P {0) - Xl3P (1) - A 2 /3P (2) + 0(X 3 ) 



in terms of the densities 



Let us first evaluate 



Pa = /4 0) + Ap a 
Ap a =\p^+\ 2 p^+0(\ 3 ). 



(7.19) 



(7.20) 



(7.21) 



We have already obtained the leading term of (if (i.e. at the MF or 0-loop order) in sec (VB) as a functions of the 
densities, i.e. for the homogeneous system considered here 



0fM [{ Pa }} ee (3Jmf [{ Pa }\ = P/hs [{ Pa }} - Pa v S a ■ 



At this point we remark that 



[{ Pa + A Pa }} = /?/«» [{ Pa }} + u a A Pa ^<Ap a |C$|A P/3 } + 0(\ 3 ) 



(7.22) 



(7.23) 



where is the two-body MF proper vertex. In the case of the SPM the expression (5.22) of takes the form 



G$(l,2) = cj&O-u)- ^(1,2), 

Pa 



(7.24) 



where c^ s is the usual direct correlation function of the HS fluid at the density p(°K Reporting the expression (7.24) 



of (7^g in eq. (7.23) and taking into account the clcctroneutrality condition at the MF level one finds 



0f<® [{ Pa + A Pa }} = [{ Pa }} + v a A Pa 



= nm 



Pa 



a Pa 



+ 0(A 3 ) 



(7.25) 



Reporting now the expression (7.17) of p a X ^ in the above equation one obtains 
[{ Pa + A Pa }} = /?/«» [{p Q }] + v a Ap a 
A 2 /? 2 



+ 



A (0)2 (Q) 



o (0) o 3 



2 \ 



Pa <l a 



>+0{\ 3 ) 



(7.26) 



where we have made use of the OZ equation for the reference HS fluid c^(0) = ft^(0)/(l + p (a) h ( ^ s (0)). The first 
contribution (7.21) of (3f can thus be recast under the form 



PaVa - (3P<® [{v a }\ = [{p a }} 
A 2 /? 2 



A 2 (0) {[p a qlY h HS (Q) + Pa q a 



p a q a 



+ 0(\ 3 ) , 



(7.27) 



where the smeared propagator A r as well as the pair correlation h H g of the HS fluid can be evaluated at the total 
density p = ^ Q p a (rather than p^) at this order in A. 

The second step in the calculation of (3f is to reexpress fiP^- 1 ' in terms of the densities p a . We note that 



\ r f «-(°) 2 



-r 2 (k) 



(7.28) 

2 , .2 „(0) 2 



which shows that j3P^' is a function of the sole squared Debye number k 2 . Defining the increment An" = k~ - / 
and performing a Taylor expansion of around «(°) 2 yields 
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A/3P (1) =~J d"k log (l + ^r 2 (k)j + ^A T (0)A« 2 + 0(X 3 ) , (7.29) 
which, after substitution of the expression of Ak 2 at order 0(A) (cf eq. (7.18b) allows us to write 
A/?P (1) = ~ J d^k logj^l + ^f 2 



- *f A 2 (0) \[p aq l} 2 h H s(0)+p aq t %|L} +G(A 3 ) , (7.30) 

where, once again, the smeared propagator A T as well as the pair correlation h HS of the HS fluid can be evaluated 
at the total density p = 'J2 a p a (rather than p^) at this order in A. 

It remains to reexpress (3P^ in terms of the densities p a , which is done readily, and to gather all the intermediate 
results. After doing this job, one finds that the free energy at the second order in the loop expansion reads as 

0f [{/>«}] = Whs [{P.}} + ljd V k (log (l + ^r 2 (fc)) - ^r 2 (k) 
-^[P°&\ 2 J d s rh HS , P (r)Al K (r) 

+ ^ W] 2 / d 3 r K (r) , (7.31) 

where a subscript p has been added to the function has, p( r ) to emphasize that it has to be computed at the density 
p = ^2 a Pa- There are good theories for the HS pair correlations fiHs( r ), f° r instance one could try hns{r) — hpy{r), 
which makes the two-loop expression for the free energy manageable. Similarly the smeared propagator in the RHS 
of eq. (7.31) must be computed from (7.13) with n 2 — 47r/3p Q <7 2 . 
Several comment are in order at this point. 

(a) We note that only the irreducible diagrams D3 and D5 of fig. 3 have survived to the Legendre transform; note 
that, moreover, the symmetry factor of D3 has changed as the result of many compensations of equal terms. Diagram 
D3 is interesting because it can be interpreted as the contribution of an effective attractive interaction — A 2 K (r) to 
the free energy. Such a term should play an important role in the description of the liquid-vapor transition of the 
SPM. 

(b) Digging the literature we have found that the expression (7.31) when specialized to the case of the RPM (i.e. 
with the last term of the right hand side set to zero) is not new and has been obtained more than thirty years ago 
by Chandler and Andersen 15 in the framework of the mode expansion formalism. According to the numerical study 
performed by these authors eq. (7.31) gives reasonable results in the liquid regime. 

(c) Recent developments of the theory of collective variables (in its modern formulation 51,52 ) yield, in the case of 
the RPM, to an expression for (3f that disagrees with eq. (7.31). As far as I understand the details of ref. 48 one 
can reconciliate the two points of view by replacing fins, p(r) by its Fourier transform at k — in eq. (7.31). This 
approximation is correct if the correlation length associated to charge fluctuations (i.e. £c ~ V K ) 1S l ar g e compared 
to the correlation length associated with the HS fluid which does not seem to be fully justified in general. 

(d) We have checked in the case of the RPM that, expanding the expression (7.31) of f3f in powers of (3 with the 
help of eq. (7.14) gives back the high-T expansion discussed in refs 3 , i.e. in reduced units (x = r/cr, j3* = /3q 2 /<r 7 and 
P* = P* 3 ) 

P\fnPM(pl=f3*f H s(pn-p*lo g 2- 2 -^( P *n 3/2 - I <P£ hHS £ {x) 

+ tt 1 / V/T) 5 / 2 J d 3 x hHS 'P'^ + (P* 3 ) . (7.32) 

This result corroborates the conclusions of Stell in his study of the relation between the 7-ordering and the mode 
expansion. 00 

VIII. CONCLUSION 

The KSSHE field theoretical representation of liquids provides a general framework for studying cither neutral 
atomic liquids 26 or ionic solutions, as in this series of papers. It is taylor made for building perturbation theories with 



27 



respect to a reference fluid chosen conveniently in general as the HS fluid. The technics is roughly always the same, 
i.e. an ordering of the cumulant expansion of the grand potential in ascending powers of some small parameter. The 
latter can be either a physical parameter such as the fugacity or the inverse temperature as considered in papers 1 
and II or, as in the present work, an abstract one related to the numbers of loops of the Feynman diagrams retained 
in the expansion. The salient features of this loop expansion can be summarized as follows. 

• The zero- loop order approximation of the free energy (3Amf constitutes a rigorous lower bound for the exact GC 
free energy (3 A. An optimized MF theory can be obtained by maximizing (3Amf with respect to the smearing 
distribution functions in the cores. The MF grand potential logS MJ ? constitutes a rigorous upper bound for the 
exact grand potential. Both functional (3Amf[{Po}] and logSMF^a}] are strictly convex in the fluid phase 
which rules out a fluid-fluid phase transition at the MF level. 

• At the MF level (or in the Gaussian approximation) the pair correlation functions coincide with those of the 
RPA theory. The direct correlation functions of the PM in the optimized MF theory are very similar to those 
considered in the ORPA (or MSA) approximation of the theory of liquids, however the pair correlation functions 
gMF, ap{r), which have simple analytical expressions, do not vanish in the cores, except in limit cases. 

• The one-loop free energy is identical with that obtained in the RPA theory of liquids. 1 ' 21 ' 37 

• An explicit and manageable expression of the two-loop order free energy can be written for the SPM; in the 
case of the RPM it coincides with a result obtained by Chandler and Andersen in the framework of the mode 
expansion theory. 15 However the expression (7. II) of the pressure seems to be a new result even for the RPM. 

The homogeneous specific free energy (3f of the PM can be used to study the critical point (CP) of the PM at 
the MF level. Recently, the one-loop order expression of fif in the WCA scheme was considered to study the CP of 
the RPM. 29 ' 48 We have performed a similar study for the SPM with the one- and two-loop order expressions of j3f 
derived in this paper. These results will be discussed elsewhere. 
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APPENDIX A: HIERARCHIES 

In this appendix we establish the hierarchies for the various distribution functions introduced in sec. (Ill A). 

a. Hierarchy for the G^ T 

As a trivial consequence of the definition (3.5) of the generator 6(1) and of eq. (3.7) we have 

Q(n + 1) T (l, . . . , n) = t (1, . . . , n + 1) , (Al) 

which is valid for all n > with the convention that G^ -0 ^ T = log 5. 

b. Hierarchy for the 
Let us apply the linear operator 8(1) to eq. (3.6). This gives 



Q(n + 1) (1, . . . , n) = */3 1/2 [G^ L > (1, . . . , n + 1) - p c (n + 1) G%> (1, . . . , n) , j (A2) 

where pc = g[? -1 ^ is the equilibrium charge density. The relations (A2) are valid for all n > with the convention 
that g£ =0) = 1. 

c. Hierarchy for the G^ 
Let us first rewrite the definition (3.10a) of more explicitly i.e. 

GW(l,...,n) = S- 1 W- 1 J V<p<p(l)...<p(n) * 

x E HS [{T7 a + i(f> a }] exp ^-1 (^v" 1 ^)) • ( A3 ) 
Then we apply 0(n + 1) to both sides of eq. (A3), which yields 

0(n + 1) GW(1, • • • , n) = -G< n »(l, • • • , n) E'^n + 1)3 [{u a }] 
+ H-W" 1 1 P<^(l)...^(n) x 

x exp ^-i (^v" 1 ^)) ©(« + 1) S HS + • (A4) 

Taking advantage of eqs. (3.6) (for n = 1), (3.21), and of the relation (3.19) one concludes that 
0(n + 1) GW(1, ...,n) = GW(l,.., n )^ A n+i < ¥>(« + 1) > n + 

+ S- 1 A/"" 1 | 2V p(l) . . . <p( n ) exp (-1 <^-V>) • (AS) 

Now we perform a functional integration by parts which gives us 

0(n + 1) GW(1, . . . ,n) = =^ [g^+^I, . . . , n + 1) - G^(l, . . . , n) < ip(n + 1) >„" 

n 

5(n + l,j)G^- 1 )(l,...,i-l,j + l,...,n) , (A6) 
which is valid for all n > 1 with the convention G^ = 1. 
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d. Hierarchy for the T 



Let us apply eq. (A6) in the case n = 1; we get 

A 2 



9(2) G«(l) = ^[< ¥>(1) >«< ¥>(2) >« -GW(1,2)J -5(1,2) , 
= ^^ a)T (l,2)-«(l,2). 



A similar brut force calculation yields 



6(3) GP> T (1,2) = -^G£> T (1,2,3). 



(A7) 



(A8) 



We want to prove that eq. (A8) is valid for all n > 3. The proof will be by induction. Suppose that, for all 2 < m < n 
with n > 3 we have either (A7) or (A8) if to = 2 or to = 3 respectively or, in other cases (m ^ 2, 3) 



6(to) G^ T (l,...,m-l) 



-A r . 



V ^ ' 47;- V 

We want to prove (A9) for m = n + 1. From the definition 



G^ T (1,...,to). 



(A9) 



(A10) 



where the sum runs over all the partitions of the set {1, . . . , n}. Now apply 6(n + 1) to both sides. One obtains 

e(n+l)GW T (l «) =1-11, (All) 

where 7 = 6(n + 1) G^ n) (1, . . . , n) is given by eq. (A6) and II results from the application of 0(n + 1) upon the sum 
of products of the RHS of eq. (A10). Since only functions T (l, . . . ,to) of order to < n are involved in these 
products, we can apply eq. (A9). Therefore 77 is made of all the terms of 



47T 



m<n+l 



47T 



G("+ 1 >(l,...,n + l)-G("+ 1 ) T (l,...,n+l) 



(A12) 



except : 

• (a) the term 



.^± <(p (n+l)> H G(") T (l,...,n) 



(A13) 



(b) the terms involving < ip(n + 1) >u, the sum of which is given by 



4-rr 
A„+i 

47T 



(m)T 



<^(n+i)> w E ro <„n^ r 



4 n) (l,...,n)-4 B)T (l,...,n) 



(A14) 



Moreover, additional terms are generated in 77 in the case where 6(n + 1) acts on some < ip{j) > (1 < j < n), this 
correspond to the special case (A7). The total contribution of these terms is given by 



f(m) T 



G 



(n-1) T 



(1,.. ., j - 1, j + 1,.. .,n) 



- E?=i <*(« + 1. j) 4 n_1) (i. ••••./ i-./ • i 

Gathering the intermediate results, one finds for 77 



(A15) 
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47T 

<<p(n+l) >n G<">(l,...,n) 



" = [G$T +1) (1, ■ • • , n + 1) - G<r +1 ) r (l, . . . , n + 1) 



- 5(n + 1, j) G£-V (1, . . . , j - 1, j + 1, . . . , „) . (A16) 

Finally we conclude that 

6(n + 1) (?(») T (l, . . . , n) = I II = ^tl G (n+i) r (1> . . . , n + 1} , (A17) 
We have just checked that if eq. (A9) is valid for m < n it is also valid for m — n + 1, therefore it is valid for all n. 

APPENDIX B: A DERIVATION OF THE FIRST SL RULE 

We derive here the first SL rule (4.9a) in the KSSHE formalism. Let us define the two-points function 

A(l, 2) = H-W" 1 J V V cxp (-1 (^"V)) ^y2 ffs [{F a + ^4] . (Bl) 

First, we integrate by parts; it gives 

. r A T / 1 . . M _ 

PffS, a (2) 



A(l,2) = -S-W- 1 / ^ [ex P (-^ (^"V)) 



= ~Wl)PffS,a(2))„. (B2) 



Another way to evaluate the function A(l, 2) is to make use of the identity (3.19) in eq. (Bl). This gives 
4(1,2) = E- 1 Af- 1 J Vy exp^-i^" 1 1^6(1) J^fns [{V a + i<t> a }] 

= l [3 1 / 2 ~- 1 M- 1 J cxp ^—(^v-^Ehs qpMl, l')GP Sifla (l', 2) 



ipV> q/} T (l,l')Gfl(l>,2), (B3) 



where we have made use of eq. (3.12). Comparing eqs. (B2) and (B3) we thus have the identity 



Moreover, we already know that 



<p(l) PHS , a (2)) n = ^ 1/2 9/3 70(1, 1') Gg(l', 2) . (B4) 



(V(l)>w = ^ 1/2 IP rp(l, 1') p /3 (l') (B5a) 

(PHS,a&)) n =P«(2) ■ (B5b) 
Taking the product of (B5a) and (B5b) and substracting the eq. (B4) yields 

Ml Ws, «(2)>£ = */? 1/2 <Z, r,(l, 1') G>g> T (l', 2) . (B6) 



We now specialize to the case of a homogeneous system and consider this equation in Fourier space in the limit k — ► 0. 
We assume that the field correlation in the LHS of eq. (B6) is regular in this limit which gives 

^Gg T (fc = 0) = 0, (B7) 

Which is equivalent, in terms of the pair distributions hp a , to the first SL rule (4.9a). 
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APPENDIX C: THE ELECTRONEUTRALITY AT THE ONE-LOOP ORDER. 



In the general case, and at the one-loop order, the densities of species a reads as 

p Q (l) = pi 0) (l) + Api 1 )(l) + O(A 2 ) 
M , SIosATa- 

'Pn—stif- (C1) 

We have shown (cf eq. (5.10)) that the MF densities p^ = Pmf. a of the homogeneous system satisfy to the elec- 
troneutrality condition. In this appendix we prove that it remains true at the one-loop order. 
We have 



= -^ (2 ' 3) ^i) A ^ 1(2 ' 3) ' (C2) 

where we have applied Wick's theorem (tp(l)(p(2)) a- = A^(l, 2). It follows from the definition (6.2) of the propagator 



A^ that 



p«(l) = -^(2,3)-^G^ T c [{F 7 + z0 7 }] (2,3) . (C3) 

(2) T 

The functional derivative of G HS c can be recast under the form 

O (2) T (3) T 

5v a (\) G HS, c( 2 ' 3 ) = G HS, ocf 1 ' 2 ' 3 ) 

- /3G^ T c (2,3,4) Uc (4,5)^r 5 (5,6)G^/ 5ct (6, 1) , (C4) 
where we have defined the hybrid density-charge correlation function 

G^ aC (l, 2, 3) = qpwp, 2')r 7 (3, 3')G$ s T af}y (l, 2', 3') . (C5) 
Combining eqs. (C3) and (C4) one finds that 

P W(1) = -^(2,3){G^ T QC (1,2,3) 

-/3G^ T c (2,3,4K(4,5)g 5 r 5 (5,6)G M F,5a(6,l)} . (C6) 
The one-loop correction to the charge density Pc\l) = q a T a (l, l')p^(l') is therefore given by 

pg>(l) = - f A^(2, 3) {G^ T C (2, 3, 1) - /3G$7 C (2, 3, 4)« c (4, 5)G^/ C (5, 1)} . (C7) 
Now it should be clear to the careful reader that 

« c (4, 5)G^ F T C (5, 1) - (V(4)pc(l)>Sc, mf , (C8) 

where V and pc are respectively the microscopic electric potential and charge density in a given configuration. With 
this remark the one-loop charge density reads as 

pgHlH-^(2,3){G^ T c (2,3,l)^ . (C9) 

In sec (V) we made the remark that the SL rules were satisfied in the MF approximation, in particular the Carnie-Chan 
rule (4.4) is satisfied and one has 

pj d(l) {V{A)pc{l)) T GC , MF = 1- (CIO) 
Integrating cq. (C9) over the volume of the system and making use of the Carnie-Chan rule (CIO) gives us 

/ d(l),g>(l) = (Cll) 

leading, for a homogeneous system, to the local charge neutrality condition p^ = 0, as expected. 
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APPENDIX D: THE FUNCTIONS T c FOR THE SPM 



The aim of this appendix is to determine the form of the functions G^g c (l, . . . , n) (cf. eq. (3.15)) for a homogeneous 
SPM. In this case all the smearing functions r Q are equal to the same t and we have 



U HS, C 



(l,...,n)=q ai ... q an r(l, l') . . . r(n, n') G^ Ql ... Q „(l', . . . , n') . 



(Dl) 



The functions G^ s c (l, . . . , n) are to be computed at the saddle point where we have shown that, for a homogeneous 
system, the electroneutrality condition 



Phs, 







(D2) 



holds. 

The HS reference system of the SPM is thus a mixture of M species of hard spheres with all the same diameter 
a . Two hard spheres of different species differ solely by their charge. However these charges do not contribute to 
the configurational energy and can thus be seen as mere internal degrees of freedom which allow to distinguish the 
various species. Denoting by z a — cxp(^ Q ) the activity of the species a we have therefore 



-HS [{z a }} = ^HS, 



(D3) 



where E.hs, o is the GC partition function of the usual HS fluid. Note that eq. (D3) is valid for both homogeneous 
and inhomogeneous systems. Of course, the density correlation functions of the mixture are also simply related with 
those of the pure fluid. 

Let us first consider the following type of n-body correlation functions 



(n)T rr s A 6 n logE HS [{z a }] 

Phs., ai ... an [{za}] (1, .-.,») = [J W 5 Zai {l)...z an {n) 



(D4) 



i=i 



It follows from eq. (D3) that 

n 

fi,. a , [{*,}] (i,...,n) = n 



i=l 



Phs, o(i) 



(") T r 1/1 



with self explanatory notations. Defining as usual the functions h^ n > T by the relations 

PHS T ai ...a n [{ Z «>] (!,■■ - ,n) = PHS, ai(l) • ■ ■ PHS, a„ (") ^....a. (1. ■ ■ ■ . «) 

we thus have simply 

h HS T ai ... an [{z*}} (1, ■ ■ ■ ,n) = [z] (1, . . . ,n) . 



(D5) 



(D6) 



(D7) 



Unfortunately the relations between the G^ s T ai a (cf. eq. (3.1)) of the mixture and those of the pure HS fluid 

are not as so simple as cq. (D7). They can be obtained by first establishing the relation between the G^ T and the 
p( n ) T which can be done as follows: we consider a simple fluid made of a single species of particles and write down 
the hierarchies for the G^ T and the T 

50^ (1 „ } = t 

oi>{n + 1) 



5p^ T {l n) = (n+1) _ _ n + + £ + (n) _ 
0^(n + 1) ii — ' 



(D8) 



Noting that G^™ -1 ) T = T = p and applying eqs. (D8) for n = 1, 2, etc one finds 
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T (l, 2) = T (l, 2) + 5(1, 2)p(l) , (D9a) 
G^ T (l, 2, 3) = T (l, 2, 3) + [5(3, 1) + 5(2, 1)] p( 2 ) T (l, 2) 

+ 5(1, 2)p™ T (l,3) + 5(1, 2)5(1, 3)p(l, 3) , (D9b) 

GW T (1,2,3,4) = ... . (D9c) 

(We have omitted to report the cumbersome expression of G^ T ). The generalization to mixtures is straightforward 

G { ^Z(l, 2) = p$£(l, 2) + 5(1, 2)5 aia2 p ai (1) , (DlOa) 

Ga}a 2 a 3 (1) 2, 3) = Paia 2 a 3 ' ^, 3) 

+ 5(l,2)5 Qia2 p£ 2 2 )^(2,3) 
+ 5(1,3)5^^(3,2) 
+ 5(2,3)5^^(3,1) 

+ 5(1, 2)5 Qia2 5(l, 3)5 QlQ3 p Ql (1) , (DlOb) 
0^3^(1,2,3,4) = .... (DIOc) 

Returning to our mixture of hard spheres one applies the eqs. (D10) by taking into account the relations (D6) and 
(D7). For instance for n = 2 this yields 

G [ H ) s T aia2 (l,2)=pHS, ai PHS., a2 hf s T Q (l,2)+p HS , ai 5(l,2)5 aia2 ■ (Dll) 

Reporting the expression (Dll) of G y H ' s a±a2 in the definition (Dl) of G y H ' s c one finds 

G^ s T c (l,2) = p HS , a q 2 a 5(l,2), (D12) 

where we have made use of the electroneutrality condition (D2) and of the normalization condition t(0) = 1. The 
same job can be done for G^ S T C and G^ S T C leading, after some tedious algebra, to the relations (7.8) reported in 
the text. 
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